New Integrable Sectors in Skyrme and 4-dimensional 

CP 71 Model 



C. Adam a )* J. Sanchez- Guillen ^ and A. Wereszczynski 6 ^ 

a )Departamento de Fisica de Particulas, Universidad de Santiago 
and Instituto Galego de Fisica de Altas Enerxias (IGFAE) 
E-15782 Santiago de Compostela, Spain 

^Institute of Physics, Jagiellonian University, 
Reymonta 4, 30-059 Krakow, Poland 

February 2, 2008 



Abstract 

The application of a weak integrability concept to the Skyrme and 
(jpn mo dels i n 4 dimensions is investigated. A new integrable sub- 
system of the Skyrme model, allowing also for non-holomorphic so- 
lutions, is derived. This procedure can be applied to the massive 
Skyrme model, as well. Moreover, an example of a family of chiral 
Lagrangians providing exact, finite energy Skyrme- like solitons with 
arbitrary value of the topological charge, is given. In the case of CP n 
models a tower of integrable subsystems is obtained. In particular, 
in (2+1) dimensions a one-to-one correspondence between the stan- 
dard integrable submodel and the BPS sector is proved. Additionally, 
it is shown that weak integrable submodels allow also for non-BPS 
solutions. Geometric as well as algebraic interpretations of the inte- 
grability conditions are also given. 
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1 Introduction 



Among nonlinear field theories, many models with topological solitons ap- 
pear to possess particular importance in various branches of physics (see for 
example pQ, [2|, [H|). In (1+1) dimensions one can take advantage of the 
well-known methods to study the mathematical structure and dynamics of 
such objects (called kinks or alternatively domain walls), like the inverse 
scattering method, Lax pair formalism or the Backhand transformation ||], 
[H]. In principle, these methods are connected with the integrability property 
of models, which implies the existence of infinitely many conserved currents. 
As a consequence, the chances for the construction of exact, topological solu- 
tions are strongly correlated with the integrability. Unfortunately in higher 
spatial dimensions there are no such general tools allowing for a systematical 
construction of solitons. What one has instead are integrable subsectors, like 
the holomorphic solitons of Belavin and Polyakov in 3 dimensions, predeces- 
sors of the selfdual instantons of Yang-Mills in 4 dimensions. 

However, there is a new promising approach, based on the construction 
of local zero curvature representations for nonlinear models, which enables 
one to calculate sufficient conditions for the existence of integrable subsys- 
tems, even though the full model does not have neccessarily such a property 
(which would be thereby exhibited if present j^j, |Zj)- Here, integrability is 
understood as possessing a zero curvature formulation and the existence of 
an infinite number of conserved currents. The significance of this method 
originates in the expectation that also in higher dimensional theories inte- 
grability is an important step on the way of deriving analytical solutions. 
That this works has been shown in investigations of especially useful models 
with Hopf solitons, where an integrable sector of the Nicole model has been 
found jS]. It was also proved that the simplest Hopf soliton belongs to this 
sector. Further, a model with infinitely many exact hopfions has been given 

i, nm. 

As is well known, a general sufficient condition (so called generalized BPS) 
for the Skyrme model |H]-PS1 analogous to the Baby Skyrme or self-dual 
cases, has not been found. Simplified models based on the Skyrme idea of 
avoiding the scaling instabilities, like the ones mentioned above, have been 
investigated. Among them, the interesting Skyrme-Faddeev-Niemi model. 
Unluckily the integrable sector of this well-known model [H], ^H] seems to 
contain no topological solitons. In order to cope with this problem, one natu- 
ral option is to relax the pertinent integrability condition (which in this case 
is the eikonal equation ^Bj, |E|)- The construction of such weaker integra- 
bility conditions has been recently described in [IBj- As one might anticipate, 
the space of solutions of this new, weaker integrable submodel is considerably 
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larger than the standard one. However, the problem whether it is populated 
by hopfions is still an open question. 

In the case of the Skyrme model the zero curvature method also allowed to 
discover a special integrable subsystem ^H]- The skyrmion with unit topolog- 
ical charge is a solution of this restricted submodel, but solitons with higher 
value of the baryon number do not belong to the integrable submodel. 
The main aim of the present paper is to define, by means of weaker integrabil- 
ity conditions, a new ('weaker') integrable sector of the Skyrme theory with 
a substantially richer set of solutions. It is carried out in Section 2, where, 
in addition, we generalize the procedure of Aratyn, Ferreira and Zimerman 
and introduce a Skyrme-like model which possesses infinitely many exact, 
chiral solitons with arbitrary topological charge. Using this toy-model one 
can study in an analytical way chiral solitons, their energies and profiles. It 
is an alternative approach to the so-called Baby Skyrme models |2Hj, [21] 
which also allows for analytical treatment of the solitons (baby skyrmions) 
but in the lower dimensional space-time. 

In Section 3, we present the investigations of the CP n model in any dimen- 
sion with similar results. In particular, in the case of the CP n model in two 
spatial dimensions, we discuss connections between strong/weak integrability 
and BPS/non-BPS sectors. The conclusions are summarized in Section 4. 



2 Skyrme model 

2.1 Standard integrable submodel 

The Skyrme model, without potential term for the chiral field, is given by 
the following formula 



L = f -fTr (u^d^Utf d^U) - -^Tr [tfd^U, tfd^U) ( 1 



A v * ) 32e 

where /„., e are constants and U is a SU(2)-vahied matrix field parameterized 
in the standard manner as 

U = e^ T \ (2) 

Here, 7$, % — 1,2,3 are the Pauli matrices and & are real fields. However, 
using results of Reference ^Hl it is convenient to take advantage of the slightly 
different parametrization 

U = e* T . (3) 
Here 



z = V£i + ei + a (4) 
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and 

^_ 1 ( \u\ 2 -l -2iu \ 



where the complex field appears due to the standard stereographic projection 

t = TTTli ( M + M *' " M 2 " X ) • ( 6 ) 
£ 1 + |up v y 

Then, in the notation of fH] we derive the equations of motion 

D»B ll = & t B ll + [A' t ,B lt ) = (7) 

where 

A,, = - ^ ^-i<9 M wr + - idpifr- + -{ud^u* - u*9 M u)r 3 

B„ = -z^r 3 + (e^r + - e^r.) , (9) 

whereas 



1 + \u\ 2 

sin 2 f 

R » = ^ - 8A (1 + | J 2)2 (iv, + iv;) (10) 

5, = M « + 4A [m, - ( ^| 2 ) 2 ^) (11) 



and 



= {d v ud p u*)d^u - (d„u)%u* (12) 
= {d v d v i)d^ - (duu)%u (13) 
= (d v ud v u*)d^ - {d v id v u)d^u*. (14) 
In addition, t± = (t\ ± zr 2 )/2. 

It has been proved JH] that the standard integrable sector of the Skyrme 
model is defined by imposing two constraints 

S^u = 0, R^u = (15) 

or in a more restricted form 

(<Vz) 2 = 0, 8^8^ = 0. (16) 

One can notice that the first integrability condition is nothing else but the 
complex eikonal equation. 

Then, one can construct two classes of infinitely many currents 



q dG dG . . 



j{H t , Hi ) =4sin ^ C0S ^ HlSfl + H 2 S;) - (1 + M 2 ) 2 + ^ J2„ (18) 

where G is an arbitrary function of £, w, m* whereas Hi, H2 depend only on u 
and it*. Indeed, such currents are conserved if the fields satisfy the integra- 
bility condition ([150. The importance of this integrable submodel becomes 
transparent if one notices that the skyrmions with charges Q = ±1 as well as 
the rational map Ansatz, which is widely used to approximate true numeric 
solitons (22]-|2H|, belong to it. 

On the other hand, it is known that the eikonal equation constrains the space 
of solutions in a rather considerable way. For example, as was mentioned 
before, there has been found only one topological soliton in the integrable 
submodel of the Nicole model. In the Faddeev-Skyrme-Niemi model the in- 
tegrable submodel seems to contain no solitons. 

In the case of the Skyrme model the assumption of the eikonal equation as 
a constraint results in a restriction of the form of the u field. Concretely, for 
the separation of variable ansatz £ = £(r) and u = u(8, <fi) a u field obeying 
the eikonal equation must be a (anti)holomorphic function of the variable 
z = tan(|)e 11 ^. If we combine it with the requirement of the finiteness of 
the topological charge, then we find that u must be just a rational map in z 
and, as a consequence, skyrmions with higher charges do not belong to this 
integrable sector. Because of that, it is reasonable to seek other integrable 
submodels, which are defined by weaker integrability conditions. 

2.2 New integrable submodel 

It has been recently demonstrated that for all nonlinear theories in Min- 
kowski space with two-dimensional target space there exists an integrability 
condition which is weaker than the eikonal equation. Here, we will apply this 
result in the context of the Skyrme model. 
We begin with the only assumption that 

dufcdPu = (19) 

whereas, contrary to the standard integrability subsector, (du) 2 is allowed to 
take arbitrary, in particular nonzero, values. Thus, we get only the second 
constraint in equation (JT5J) 

R^u = 0. (20) 
The pertinent equations of motion read 

= —^—u*S»d,u (21) 
1 + w 
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and 

= 4^j^5^X- (22) 

Since the eikonal equation is not imposed, S^d^u ^ and the left hand side 
in (|21|) is not identically zero. One can find that 

S^u = (du) 2 (l - 4A 2 (90 2 ) (23) 

The first class of conserved currents can be taken as follows 

Jfj, — G U S^ — G M *S^. (24) 

In spite of the fact that the form of the currents is identical to (fTTjl there is, 
however, a subtle modification. Now, an arbitrary function G is assumed to 
depend on the square of the modulus of u, i.e., G = G(uu*,t;). 
Then, the divergence reads 

d% = G uu d^uS ^ + G uu *d»u*S» + G M <9% 

- G u * u *avs; - G u , u d»us; - g^s; + G u d^s^ - G u *d^s; (25) 

Taking into account ()2()j) and observing that 

s»d^ = 0, dvs,, = d^us; 

we get 

d^j^ = G"{u* 2 d^uS^ - u 2 d^ i u*s;) + G'{u*d^s^ - ud^s;), (26) 

where prime denotes differentiation with respect to the modulus squared. 
Using (|2Tj) and (}2"3*|) one can rewrite it as 

,9% = G"(l - AX 2 (d0 2 )[u* 2 (du) 2 - u 2 (du*) 2 } + 

G' - + 2 |w|2 (1 - 4A 2 (<9£) 2 ) [u* 2 (du) 2 - u 2 (du*) 2 ]. (27) 
Thus, the current is conserved if we assume the weak integrability condition 

[u* 2 (du) 2 - u 2 (du*) 2 } = 0. (28) 
Similarly, one can consider the second type of currents, i.e., 

Jf = 4 sine cos + H 2 S;) - (1 + \u\ 2 ) 2 (|| + ^) R p . (29) 



It can be easily shown that its divergency vanishes if the functions Hi and 
H 2 are of the form 

_dh dh 
n\ - n 2 - , 

ou ou* 

where h is any function of the modulus squared uu*. It strongly simplifies 
the currents. Namely, they are reduced to the first class currents J G . 
Let us now summarize the results obtained so far. We have defined a new 
integrable sector of the Skyrme model which consists of two integrability 
constraints 

d ll £d*u = 0, u*\duf -u 2 (du*f = (30) 

and two dynamical equations (j21jl . (1221) . The pertinent infinite family of 
conserved currents is parameterized by any function G, which depends on 
uu* and £ in an arbitrary manner 

J p = G U S^ — G U *S^. (31) 

The geometric meaning of the integrability conditions ()30|) becomes especially 
well visible if we express the complex field in terms of the scalars S, A 

u = e s+iA . (32) 

Then, using equations (jSOjl we get 

8^0^ = 0, d„Ad^ = 0, d /i 'E& t £ = 0. (33) 

In other words, in the integrable sector all gradients of the scalar fields must 
be mutually perpendicular. 

It is worth to notice that such an integrable submodel can be constructed 
for a massive generalization of the Skyrme theory |2Sj as well. For example, 
introducing the most typical massive term 

V = f 2 m 2 Tr(l - U) (34) 

results in a modification of equation 



= ^ T^iZi ^ 9 ^* + 2m 2 sin^ (35) 
(1 + \u\ z ) z 

One can check that it does not influence the definition of the integrable 
sector and the conserved currents. This may be of some interest since the 
Skyrme model with a massive term serves considerably better as a model of 
the low energy QCD, describing low energy degrees of freedom i.e. hadrons 
and nucleons with acceptable accuracy j2Z|, |2~%] . 
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As was discussed in the previous subsection, all rational maps (or generally 
holomorphic functions) obey the standard integrability condition. Obviously, 
they satisfy the weaker condition as well. However, this weaker condition 
allows for a much larger class of solutions. For instance, non-holomorphic 
functions in the form ()32j) fulfill the constraint. It may be of some interest 
since it has been recently shown that a non-holomorphic Ansatz approximates 
true Skyrme solitons considerably better than the standard rational Ansatz 
|29j . [80] . In particular, the non-holomorphic approximation of the skyrmion 
with topological charge Q = 2 [29J obeys our integrability condition. 



2.3 Integrable chiral model 

One might ask whether it is possible to further relax the integrability condi- 
tions (J3()j) . In fact, one can construct a chiral SU{2) model which is integrable 
even if the second constraint in ()30|) is neglected. In other words, there are 
no additional requirements for the the complex field u. This property, as it 
will be shown below, allows for the existence of solitons with an arbitrary 
value of the topological charge. 

The model we are going to focus on is given by the following Lagrangian 
density 

L = -m)9(u,u*) [K^uf* + [(cU) 2 ] § , (36) 

where / and g are arbitrary, differentiable functions depending on £ and u, u* 
respectively. Additionally 

K„ = (d v ud u u*)d^u - (d v ud u u)d^u*. (37) 

The fractional exponent presented in the upper formula is understood as 
follows 

[K^n*]* = K^u*\K^u*\-^. (38) 

It can be observed that such a form of the model guarantees also an inter- 
esting way of circumventing Derrick's theorem since the energy is invariant 
under scaling transformations. This attempt was originally proposed 30 year 
ago by Deser et al 31 to consider the theory of pion fields n l and further 
developed by many authors (see Nicole model [H], Aratyn-Ferreira-Zimerman 
model jinj and their generalizations jH2], jSH])- Let us notice that our model 
can be treated as the modified AFZ model coupled in a non-minimal way 
(via the dielectric function /) with a non-standard scalar field 1 . 

1 Models with a non-canonical kinetic term for the scalar field have been recently widely 
discussed in the cosmology. They have been applied in the context of inflation 34. , dark 
matter [35] as well as the modified Newtonian dynamics (MOND) 36 . For the detailed 
studies of the kinetic term occurring in our model (|36[) see |37| . 
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In principle, one may rewrite this model in terms of the original chiral field 
U but in practice such a Lagrangian would take a completely illegible form. 
The pertinent equation of motion for the complex scalar field reads 



3 x 
2 M 



fg [K^uT 1 K»\ - g u *f [K^u*]* = (39) 
where g u * = d u *g. This equation can be simplified to the following expression 



d„ 



fg* [K^uT' 4 K» 







or in the most compact form 



where 



= 0, 



(40) 



(41) 



= fg* [K^u*]" K*. (42) 
The second independent equation of motion, for the real scalar field £, is 



3d,, 



(43) 



In order to establish the integrability property of the model we introduce the 
following currents 

Jix — fc/j,G u — fc*^G u *, (44) 

where G is an arbitrary function of u, u* and £. Taking into account the 
equations of motion it can be proved that these currents are conserved if we 
assume only one constraint (fl9"jh It is due to the fact that the expression 
JC^d^u = is a mathematical identity which does not restrict the form of 
the complex field. 

In order to find topologically nontrivial solutions of these dynamical equa- 
tions we will consider only time-independent configurations and assume the 
Ansatz for the fields 



u = u 



(</>, 9) = e in4, h(9) 



(45) 



and 



£ = £(r), (46) 

where the spherical coordinates (r, 9, <p) have been introduced. Here n is an 
integer number. Of course, such an Ansatz obeys the integrability condition. 
Thus, the obtained solitons will belong to the integrable submodel. 
One can easily find that 



K ■ Vu* = 4 



nhliQ 
sin 9r 2 



(47) 
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and as a consequence 



ic = V29 1 / 3 / 



nhho 
sin 9 



2 e im 



nh 

0, -r-z,lh 6 
sm9 



(46 



Substituting these formulas into the pertinent field equation we get the fol- 
lowing second order, ordinary differential equation for the function h 



nhg 



1/3 



nhhg 
sin 9 



nh e g 1/3 (h) 



nhhg 
sin# 



0. 



(49) 



It can be simplified to this expression 

hhf, 



do 



9 



1/3 



sin 6 



(50) 



which possesses the obvious solution 

hhf, 



9 



1/3 



sin 9 



/i. 



(51) 



where \i is a positive constant. 

Taking into account the Ansatz and the upper obtained solution (JoTj) we are 
able to rewrite the equation of motion for £ in the form 



-d r 



(2n) 3 / 2 /i 3 



0. 



After introducing a new variable 



we derive the following equation 



x = lnr 



Fortunately it can be integrated for any function /. The solution is 

(2n) 3 / 2 



(52) 



(53) 



(54) 



(55) 



Both formulas (JSTj) as well as (jB3|) can be integrated, at least using stan- 
dard numerical methods, for all reasonable functions g and / leading to the 
following general solutions 



J g 2 ^hdh = —jj 2 cos 9 + a 



(56) 
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and 

lm = ^r^ (57) 



v"7 ^2 ' /•„• 
where ao,r are integration constants. 

Let us now analyze in detail a particular case with the following forms of the 
coupling functions 

9( "-"* ) = (TTirpF (58) 

and 

/(0 = sin 3 £. (59) 
The corresponding Lagrangian density can be expressed as follows 

L — — sin 3 i [n ■ (d^n x d v n)\ § + [(<9 M 2 ] 1 , (60) 

— # 

where n = £/£ is a three component unit vector. 
Then 

(TTW^ (61) 

Hence, the solution is given as follows 

h 2 = (1 + a + /2 2 ) sin 2 | + (1 + a - /i 2 ) cos 2 f 
(— ceo — Ai 2 ) sin 2 | + (— «o + /i 2 ) cos 2 f 

This solution reduces to a well-known form if the parameters satisfied 

1 + a - ^ 2 = fi 2 + a = 0. (63) 

That is 

«o = ~ /i 2 = ^- (64) 

Indeed, then we get 

ft, = tan -. (65) 
In other words we obtained the following complex scalar field 

u{9,(f>) = e in *tan^. (66) 

One can observe that if we take into account the stereographic projection 
then we find 

n = (cos n(f> sin 9, sin n(f> sin 9, cos 9) . (67) 
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The remaining field equation 



n 



^2 



sin£ 



(6* 



can be solved as well, providing the solution 

T \ ^" 

£ = 2 arctan ( — 

^0- 



(69) 



The topological charge of the obtained configuration can be computed from 
the standard expression 



Q 



1 f ^.du A du* A d£ 



An 2 J (1 + M 2 ) 2 ' 
Inserting our solutions (JfifiJ) . (JfiHj) into (J70*|) one gets that 



(70) 



Q = nJ 



dt e(o)-e(oo) 

— = n 

7T 7T 



-n. 



(71) 



Therefore the solution represents a chiral field with arbitrary value of the 
baryon number. 

Let us now compute the corresponding total energy. In the case of static 
configurations it is given by the formula 



E 



KVu* 



+ 



d 3 x. 



In particular, for the example considered above, it can be rewritten as 



E 



n\ 2 sin 3 f 

^3T + (« 



;> 3 



r 2 sin 9d(f>d9dr. 



(72) 



(73) 



Taking into account equation (joTj) we get 



E = 4vr • 3 / drr z r. 
Jo 



Finally, using solution one finds that 

p 37r2 



(74) 



(75) 



As one could expect the energy of the obtained configurations is proportional 
to the topological charge. 
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Remark: It is easy to observe that solutions of the model (jBTIj) with other 
coupling function g(u,u*) have the same topological properties. Of course, 
it is true only if the pertinent boundary conditions for h and £ are fulfilled. 
Namely, 

f or 7r ifr^O , ^ j or oo if ^ — ^ ,„ . 

) 7r or if r — > oo ' [ooorO if # ^ oo . 

For instance, if one considers the function g in the following form 

9M = (i + \u\*y (77) 

where p > 3/2. Then the general solution reads 

3 

/ 3 \ 2 p- 3 

' / " ^ -1. (78) 



\2(2p- 3) (fi 2 cos8 - a ) J 
The topologically interesting configuration is obtained if the constants 

3-2p 

^ = 4(2^3)' ao = -^ (79) 

Indeed, 

V cos + 1/ 

is a function which satisfies the required boundary conditions (J75j) . The 
energy-charge relation remains unchanged. Assuming other forms of the 
function g it is possible to construct more complicated but still topologically 
nontrivial solutions. 

Remark: In order to derive solutions with fractional topological charge one 
can consider the following coupling function 

/(0 = sin 3 (^j , (81) 

where q is a positive parameter. The function g is assumed in the standard 
form (JH%|) providing the previously described complex field (jfifijl . However, 
the real field £ is not longer a map onto the segment [0, tt]. Indeed, now it is 
given by the expression 

£ = 2garctan (^) g ^ (82) 
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i.e., £ G [0, <77r]. Therefore, the original chiral field U does not cover the 
whole target space providing an arbitrary, in general noninteger, value of the 
pertinent topological index 

Q = -nq. (83) 
Nonetheless, the energy is still finite 

E = 3</2ir 2 nq. (84) 

As one sees, the energy is proportional to Q even if it takes noninteger value. 
The existence of finite energy solutions with any noninteger value of the 
topological index can indicate the the obtained chiral solitons are unstable. 
This is in accordance with Ref. [SZj where we have shown that the static 
Hopf solitons in the Aratyn-Ferreira-Zimerman model have analogous stabil- 
ity problems. 



3 CP n model 

3.1 Standard integrable submodel 

Let us now investigate the CP n model in an arbitrary dimensional space-time 

0, EE], E29, M 

(1 + v) ■ u)(dnU^ ■ d^u) — (u^ ■ <9„-u)(d'Vf ■ u) , . 

L = 77—1 ' 85 

[1+U~ • U) z 

where u is a column of n complex fields and u ' • u = u*u a , a = l..n. The 
pertinent equations of motions read 

(1 + v) ■ u)d 2 u a = 2{v) ■ d^u)d^u a . (86) 

As it was established in [Zj this model also possesses an integrable sector, 
given by the dynamical equations 

d 2 u a = 0, a = l..n (87) 

and the set of integrability conditions 

d^Uad^Ub = 0, 1 < a < b < n. (88) 

In fact, if we assume (|87|). (|5Hj) the following currents are conserved 

J a = -x—d,,Ui — Tr~-d u ut, (89) 
dui " du* 11 
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where T is an arbitrary function of it*. Applying the decomposi- 

tion of the complex field used before 

Ua = e^ a+iAa , a = l...n (90) 

we can find a geometric interpretation of the integrability conditions 

(d,A a ) 2 = (^S fe ) 2 , (91) 

d^Aad^b = 0, 1 < a, b < n. (92) 

Interestingly enough, in (2 + 1) dimensional space-time, the standard inte- 
grable submodel constitutes the BPS sector of the CP n model. That is, the 
solutions of the integrable subsystem (jHZj), (f88|) are nothing else but the very 
well-known BPS solutions of the CP n model 

U = (93) 

where / is a vector of f a = f a { z ) (anti)holomorphic functions of z = x + iy. 
Indeed, due to the conformal symmetry the general solution of the integrable 
submodel is just a collection of (anti)holomorphic functions 

u a = u a {z) or u a = u a {z). (94) 

Then, using the local U(l) invariance one can express the BPS solutions of 
the CP n model in the following form 



U = . . (95) 

VI + u a *u a V u J 

Of course, as one usually wants to deal with objects carrying finite topologi- 
cal charge the holomorphic function must be restricted to arbitrary rational 
functions. 

In more than two spatial dimensions the conformal symmetry does not occur 
and the chance for the existence of an integrable submodel with a nontrivial 
set of solutions becomes limited. It is easy to see that a general solution of 
the integrability conditions is u a = A a f(u), where A a are complex constants 
and / is any function of u, which satisfies the (d + 1) dimensional scalar 
eikonal equation 

dpUdPu* = 0. 

If one combines this result with the fact that such configurations should obey 
the dynamical equation (|87|) as well, then one may expect that the chances 
for nontrivial solutions of the integrable sector are highly reduced. It, of 
course, restricts in a quite considerable way the range of applications of the 
obtained integrable sectors. 
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3.2 New integrable submodels 

In this section we construct new integrable sectors in the CP n model. In 
general, such alternative submodels are found if much less restrictive con- 
straints are imposed. 
We consider currents similar to 



but the function T is assumed to depend on the field variables in a specific 
manner. Namely, 

T = T{u\u x + ... + u* k u k , u* k+1 u k+1 , u* n u n ). (97) 

Of course, we expect that the currents are conserved i.e. 

d»jf = 0, i = l...n. (98) 



Thus, taking into account the field equations (Joojl . we get these integrability 
conditions 

k 

^(ululd^Uad^Ub - UaUbd^uld^ul) = 0, a = l...n (99) 

6=1 
k 

Y^{K u bd^u a d^u* b - Uauld^uld^Ub) = 0, a = l...n (100) 

6=1 

u^uld^Uad^Ub — u a Ubd^u* a d^u* h = 0, b = k + 1, n a = l...n (101) 
u^Ubd^Uad^ul — UaUld^u^d^Ub = 0, b — k + 1, n a = l...n. (102) 
Once again, it is convenient to rewrite them using (|90p. Then, 



d^Aa (£^) = °> ^E a (^^ A& ) = °' a = 



(103) 



d ll A a & l E b = Q, d li 'E a & i A b = 0, a = l...n, 6 = k + 1, n. (104) 

As it was described in the Skyrme model in the previous section (or in 
the Faddeev-Niemi model in (HI), in an integrable sector gradients of the 
pertinent real scalar fields (real coordinated on the target space) must be 
perpendicular. However, it does not mean that all gradients are mutually 
perpendicular. Varying the form of the T function we get that for some of 
the gradients it is sufficient to be perpendicular to a sum of others. 
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To obtain even weaker integrability conditions, and in a consequence a sub- 
model with a larger set of solutions, one can take under consideration the 
following current 

n ( dT dT \ 

^ = Et^-^X • ( 105 ) 



i= i \ du i du * 



In this case the integrability conditions read 



E<vO (e^W) =o, (106) 



\a=l / \b=l 



a = fc + l,...,n (107) 



<9 M A a <9 M £ b + d^ad^Ab = 0, a = A; + 1, n 6 = A: + 1, n. (108) 

To conclude, there is a family of integrable sectors of CP n model. Integra- 
bility conditions vary from the strongest ones (the most restrictive) 

d^A a d^ b = 0, a, b = l...n, (109) 

when T = T(uiul, ...,«„«*) to the weakest ones (the lest restrictive) 

(e<W) (e^e*) =o, (no) 

if JF = JF^u^ + ... + u n ul) and take all possible intermediate cases. Obvi- 
ously, the full integrable submodel consists of the field equations (JHH)) and a 
choice of an integrability condition, as discussed above. 
It should be underlined that the construction of such integrable sectors of the 
CP n model remains unchanged if we add a potential term to the Lagrangian 



Let us again interpret the derived integrable sectors in the context of CP n 
model in (2 + 1) dimensions. One can check that the most restrictive con- 
straint ()109j) possesses, in addition to standard ones (JMJ), the following solu- 
tion 



A a = A a ^-j, E a = Z a {zz*), (111) 

where A a , E a are arbitrary functions depending on phase and modulus respec- 
tively. On the other hand it is known that the CP n model has also non-BPS 
solutions, which can be constructed from the BPS solitons by acting with a 
projective operator jH]-jllj 

pk f 

u= \nry (112) 
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where 

P+f = d g f-(fd z f)-£- 2 , k = 0,..,n. 

These non-BPS solutions (or at least some of them) obey the integrability 
conditions. As an example let us consider / = (1, z, z 2 ), then for k = 1 



U 



l + 4|z| 2 + 6|z| 4 + 5|z| 6 + \z\ 





z*(l + 2 


A 2 ) \ 




M 4 - 




V 


-z(2 + 


A 2 ) ) 



(113) 



Now, after expressing it by means of the parameterization (|95|) we get that 
the pertinent scalar complex functions are in the form (jlllj) . Since all non- 
BPS solutions by construction fulfill the dynamical equation of motions (j8fij) . 
we can conclude that new integrable sectors of the CP n model in two spatial 
dimensions consist of BPS as well as non-BPS states. 

As we see, in (2 + 1) dimensions, there is a striking correspondence between 
integrability and BPS property of solutions. The strong integrable sector is 
just the BSP sector whereas weaker integrable sectors contain also non BPS 
solutions. In higher dimensions we usually do not have any BPS solitons, 
however, the strong as well as weaker integrable sectors are still well-defined. 
Thus, integrability is an important tool in investigating nonlinear models in 
higher dimensions. 

An algebraic meaning of the family of integrable sector can be easily found if 
one considers volume preserving diffeomorphisms of the volume 2n-form on 
the CP 11 target space 

VL = g(uul, u n u* n )du\ A du\ A ... A du n A du* n . (H4) 

In other words we are interested in transformations u a — > v a (ui, «*, u n , m*) 
which leave the 2n-from unchanged 

Q = g(uu{, u n u* n )dui A du\ A ... A du n A du* n = 

g(vvl, v n v^)dui A dv\ A ... A dv n A dv*. (H5) 
Let us discuss infinitesimal transformation u a — > u a + e a . Then 

du a — > du a + e a Uh du b + e a u *du* b , (116) 

where summation over indices b is performed. Now, the invariance condition 
()115j) gives 

qU 1 + e *\ + ... + e n + e™) + 
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{d ig ){ Ul e* 1 + ^e 1 ) + ... + (d n g)(u n e* n + u* n e n ) = 0. (117) 
Here e° a = d Ua e a and dig = d UlU *g, etc. After introducing 

t a = g- l F K (us) 

we can rewrite equation ()117|) in a very compact form 

(d ui d ul + ... + d Un d K ) (F + F*) = 0. (119) 
The general solution reads 

F + F* = EC«%-,< 1 ), (120) 

i 

where <7j = 1,2 and u 1 = u, u 2 = u* . However, from our point of view it is 
sufficient to investigate only pure imaginary solutions 

F + F* = 0. (121) 

Thus, F = %G where G = G(ui, u\, u n , w*) is an arbitrary real function. 
It follows that all volume preserving diffeomorphisms are generated by the 
vector fields 

n 

v G = ig- 1 E (G<9u a - G Ua d K ) , (122) 

a=l 

which satisfy the Lie algebra 

n 

[v G \v G >] =v G \ G 3 = E {Gi K G 2Ua - G 1Ua G 2< ) . (123) 

a=l 

In order to find Abelian subalgebras we assume that G\^ take the following 
form 

G = G{u x u\ + ... + u k u* k , u k +iu* k+1 , u n u* n ), (124) 

where k = l..n is a fixed number. Indeed, the corresponding commutator 
vanishes. Thus for each value of k we get an Abelian subalgebra Tik- It is 
straightforward to see that these subalgebras form a tower of algebras 

U x D H 2 D ... D H n - V (125) 

Further, the vector fields v G can be identified with the Noether charges of 
the currents ()105j) . once the identification G — > JF is made. 
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4 Conclusions 



We have shown that the method, developed in [TH] and originally devoted 
to the construction of integrable submodels for S 2 valued field theories, can 
be easily adapted for models with different target space, and defined in an 
arbitrary dimensional space-time. 

The interesting point is that the new derived integrable sectors of Skyrme 
as well as CP n models are defined by integrable conditions which are con- 
siderably weaker than the standard ones. Thus, one might expect that the 
set of solutions of the integrable subsystems is larger. In particular, it may 
cure the problem of the apparent non-existence of soliton solutions with a 
topological number largen than one in the standard integrable sector of the 
Skyrme model. However, the problem whether the new integrable submodel 
possesses such topological solutions is still an open question. 
Moreover, a correspondence between BPS sector and strong integrability in 
the (2 + 1) dimensional CP n model has been established. In addition, we 
have shown that weak integrable sectors support non-BPS solutions. Such a 
relation holds only in three dimensional space-time but it may indicate that 
in higher dimensions integrable submodels can play the role of BPS sectors. 
Therefore, one can hope that the construction of integrable submodels might 
compensate for the non-existence of BPS sectors and provide a useful tool 
for constructing solitons. 
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